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Abstract A model of vesicle electrodeformation is described which obtains a 
parametrized vesicle shape by minimizing the sum of the membrane bending 
energy and the energy due to the electric field. Both the vesicle membrane and 
the aqueous media inside and outside the vesicle are treated as leaky dielectrics, 
Oh! and the vesicle itself is modelled as a nearly spherical shape enclosed within a 

thin membrane. It is demonstrated (a) that the model achieves a good quanti- 
tative agreement with the experimentally determined prolate-to-oblate transition 
frequencies in the kHz range, and (b) that the model can explain a phase dia- 
gram of shapes of giant phospholipid vesicles with respect to two parameters: the 
frequency of the applied AC electric field and the ratio of the electrical conduc- 
tivities of the aqueous media inside and outside the vesicle, explored in a recent 
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paper (S. Aranda et al, Biophys. J. 95:L19-L21, 2008). A possible use of the 
frequency-dependent shape transitions of phospholipid vesicles in conductometry 
of microliter samples is discussed. 
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1 Introduction 

An increasing awareness of the effects of the electromagnetic field on biological 
samples [ 1|, as well as its use in biotechnology [2| has instigated numerous studies 
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Fig. 1 Morphological phase diagram of lipid vesicle shapes subjected to an AC electric field. 
The two parameters spanning the diagram are the frequency of the applied AC electric field 
and the ratio of the electrical conductivities of the aqueous media inside and outside the vesicle 
(from ref. |6 |; reproduced by permission of The Royal Society of Chemistry). 



on the effects of electric field on individual biological cells. When exposed to an 
AC electric field, biological cells exhibit several connected phenomena such as 
orientation, dielectrophoresis, electrorotation and deformation J3]|4]|5l. Of these, 
deformation has probably received the least attention, despite the fact that it can 
be conveniently studied on giant unilamellar phospholipid vesicles (GUVs) J6][7), 
which serve as models of biological cells. 

Extensive studies of the behaviour of biological cells in the electric field have 
been conducted by Schwan since the 1950's [8]. However, these early studies 
treated cells as rigid objects. A theory of lipid bilayer elasticity [9| was needed 
for the first model of vesicle electrodeformation 1101 . This first model predicted 
a prolate spheroidal deformation; this prediction was later confirmed both by im- 
proved theoretical treatments 01111121 and by experiments 1131 . conducted in a 
2 kHz AC electric field. However, experiments where the frequency of the field 
was varied 11411151 have demonstrated that the vesicle shape undergoes a transi- 
tion from prolate shapes at low frequencies to oblate shapes at higher frequencies. 
Independently, such behaviour was also predicted by Hyuga and co-workers 0161 
[T71 . who described a prolate-to-oblate shape transition in certain conditions with 
respect to the electrical properties of the aqueous medium inside and outside the 
vesicle. 

A recent paper [18] extended the studies of vesicle shape in electric field to 
higher frequencies and explored the effect of varying electrical conductivity of the 
internal and the external aqueous medium in a systematic manner, thus obtaining 
a more complex morphological phase diagram (Fig.|TJ. The authors have demon- 
strated that in the case when the conductivity of the external aqueous medium 
exceeds the conductivity of the aqueous medium inside the vesicle, another shape 
transitions exist in the MHz range in addition to the previously observed prolate- 
to-oblate transition in the kHz range. 
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The rest of this paper is structured as follows: Section|2]presents the theoretical 
model, which treats the vesicle as a thin shell made of leaky dielectric and derives 
its shape by minimizing its total energy, consisting of membrane bending energy 
and the energy due to the electric field. Section [3] describes the experiments in- 
vestigating the prolate-to-oblate morphological shape transition in the kHz range. 
The results - a comparison of the model predictions with the experimental data 
for the prolate-to-oblate transition on the one hand and with the experimental mor- 
phological phase diagram (Fig. [T} on the other hand - are presented in Section [4] 
Section [5] discusses the limitations of the presented model, compares it with the 
existing ones, and introduces a possible application of frequency-dependent tran- 
sitions of vesicle shapes for conductometry. Finally, Section [6] presents the main 
conclusions. 



2 Theoretical analysis 

In general, the model follows the approach introduced by Winterhalter and Hel- 
frich [12], i.e., a parametrized vesicle shape is obtained by minimizing its total 
free energy. Two terms enter the free energy: the membrane bending energy and 
the energy due to the electric field. A closed-form expression for the membrane 
bending energy is well-known [9 |, and the change of the free energy due to the 
electric field is calculated as the work done by the force of the electric field while 
deforming the vesicle 11121 . 

Once an electric field E is applied, it introduces a single distinct axis into the 
system, so the treatment can be limited to axially symmetric shapes. This elim- 
inates the dependence on the longitude angle <j>, if the polar (z) axis is chosen 
parallel to the applied field. Thus, the vesicle surface can be parametrized as 

r(0) =s +s(9), (1) 

where \s(6) \ <C sq, sq denoting the deformation independent of the polar angle 
9. In an absence of deformation = 0), so equals the radius of the unde- 

formed sphere (ro). The deformation is independent of the sign of the electric 
field and thus proportional to E 2 in the lowest order. The field itself being pro- 
portional to cos 0, the deformation coupled with this field can be expected to be 
proportional to cos 2 6. In terms of expansion into spherical harmonics, this limits 
us to a sum of even terms. Retaining only the terms proportional to E 2 or lower, 
a quadrupolar term remains, where s(6) equals the second Legendre polynomial: 
s(9) = jS2(3cos 2 — 1), with S2 being a measure for the extent of deformation. 
Positive values of S2 indicate prolate deformation, while negative values indicate 
oblate deformation. 



2.1 Membrane bending energy 

The requirement for a local area conservation, which assures that the membrane 
stretching is independent of polar angle 9, implies that a quadrupolar displacement 
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Sr r in a radial direction is accompanied by a tangential displacement 8rg lfT2l : 

Sr r = i(3cos 2 0-l)i 2 (2) 
5re = — cos sin 0^2 • (3) 

The total membrane area expansion is determined by the relationship between sq 
and tq. Taking into account the requirement for a constant vesicle volume, the 
following relationship for sq ls obtained: 

s 2 

so = r - j 1 - ■ (4) 

The correction for a constant volume Q contains a higher term in the powers of s 2 
and thus does not affect in the lowest term either the bending energy or the energy 
due to the electric field, yielding sq = >"o an adequate approximation. 
The expression for the vesicle bending energy |9ll : 



Gbend = -jk c j> (ci + c 2 -c Q ) 2 dA 



(5) 



Here, k c is the bending elastic modulus of the membrane, while c\ and c 2 are 
the principal curvatures of the membrane. The spontaneous curvature cq vanishes 
for a bilayer composed of two equal layers. The integration is conducted over the 
total membrane area A of a quadrupolarly deformed vesicle. Up to quadratic order 
terms in s 2 , the total bending energy of a nearly spherical vesicle can be expressed 
as 10: 

487T fs 2 \ 2 

Gfeend = 87Tfc c +— k c [y\ . (6) 

This can be readily interpreted as the bending energy of a sphere plus an addition 
due to the quadrupolar deformation. 



2.2 Contribution of the electric field 

The first step in evaluating the contribution of the electric field to the total free 
energy is calculating the electric field. The present treatment is limited to small 
deviations of vesicle shape from the sphere, which greatly simplifies the treat- 
ment: instead of computing the electric field in the presence of the actual vesicle 
shape, one can compute the electric field in the presence of a spherical shell. Both 
the aqueous solution inside and outside the vesicle and the vesicle membrane are 
treated as lossy dielectrics. The aqueous solutions inside and outside the vesicle 
usually have identical dielectric permittivity, while their electrical conductivities 
can differ (Fig.[2jl. 

In order to compute the forces exerted by the electric field on the vesicle, one 
first has to compute the electric field around a vesicle. Gauss' law V • D = to- 
gether with the requirement for an irrotational electric field V x E = 0, which 
stems from Faraday's law for electromagnetic induction, yields the Laplace equa- 
tion for the electric potential U lfl9l : 



V 2 U =0 



(7) 
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Fig. 2 A vesicle is modelled in spherical coordinates (r,6) as a spherical shell with a radius r 
exposed to an external electric field E. The conductivity and permittivity of the aqueous solution 
outside are denoted by rj ex and e ex , and the corresponding quantities in the vesicle interior by 
rji n and £; n . Membrane surface transconductance and surface capacitance are denoted by g m and 
c m , respectively. 



The boundary conditions on the membrane (r = ro) require that the total sur- 
face charge density, including both free charge and displacement charge, must 
vanish. The case where a finite potential can be supported across the membrane is 
known as the series admittance limit CO, pp. 230-232). In the spherical geometry, 
this yields a system of equations 

((Tex-iffl£ex)£T Oo) = (a in - ime in )E^ 2) (r ) , (8) 
(ae X -ift)eex)^ 1) (ro) = tem-ift)c m )(fj( 2 )(ro)-[/ (1) (ro)) . (9) 

Apart from © and (0, the system is constrained by two additional conditions: one 
requiring that the electric field far away from the vesicle is unperturbed, and the 
other requiring that the electric field is finite inside the vesicle. Membrane surface 
transconductance g m and surface capacitance g m are related to membrane electric 
conductivity and dielectric permittivity: g m = O m /h and c m = £ m /h, h being the 
membrane thickness. 

In spherical coordinates, Q readily decouples into the radial and the angular 
part, and can be solved with the usual ansatz: 



£/(*) = I 

2 



lcos0e- iffl '+C.C. 



(10) 



The coefficients , b^ ;k= 1,2, which can in general be complex to allow for a 
phase shift, are determined from the boundary conditions ( 1819b . 

The conditions for an unperturbed field far away from the vesicle and a finite 
field inside the vesicle immediately yield two coefficients: 

a (1) = -£o, (11) 
b {2) =0. (12) 
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The two remaining coefficients are obtained by solving ( 1819 



7 (2) 



3E r 



Eor* 



(g m - iac m ) (a ex - itoe e , 

OexOin +gmK°ex ~ Ofo) ~ (£ex£in + C m r(e ex - £ in )) ft) 2 



- i (gm>"( e ex - Sin) + C m r(o ex - CT in ) + £i n C7ex + £ e xOi n ) ft) 



(13) 



(14) 



4 = 2a sx a in +g m r(2a ex + a in ) - (2£ ex £in + c m r(2e ex + e in )) ft) 2 

(2a ex + (Tin) + 2(£inC7ex + £exOi n )) » • (15) 

The surface density of the forces exerted on the boundary of dielectrics by the 
electric field can be computed as a scalar product of the Maxwell stress tensor and 
a vector normal to the membrane: 



f= (rW - r( 2 )) ( 



(16) 



The force vanishes in a homogeneous medium, but can in general be non-zero 
on the boundaries of media with different electrical properties. The Maxwell stress 
tensor is defined as 



T =D®E-i(D-E)7, 



(17) 



where I denotes the identity matrix. 

Unlike in the case of the bending energy term (5), for which a closed-form 
expression was obtained, an approach where energy difference is computed is 
employed here. <5Gfi e id denotes a small change in the energy due to the electric 
field, when a sphere (S2 = 0) is perturbed by a small deformation change 8s2- This 
energy difference is calculated as the work done by the forces of the electric field 
during the displacement of membrane elements <5r (equations l2l3b , integrated over 
the entire membrane area fl2l: 



5G fie id = -^(f-<5r)dA. (18) 

The integration is conducted over a sphere, which is consistent with the limit of 
small deformations (S2 <C ro)- 

Substituting the coefficients dTTHTSb into dT8l yields a lengthy expression for 
<5Gfieidi which can be written as a sum of two dispersion terms: 

5G field = -f e wE y 0Ut s 2 (§. + ^ + r JL i ) . (19) 

Three dimensionless coefficients %\, and £2, and two characteristic times %\ 
and T2 are rather lengthy expressions involving six different material constants: 

£ex> £in? Cexs ^in? Cm> gm-> 

and the vesicle radius ro. 
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2.3 Vesicle deformation 

For small deformations, quadrupolar deformation only induces small perturbative 
changes in the vesicle bending energy and the energy due to the electric field. 



G be nd(s 2 ) * G bend (s2 = 0) + \ ^gfd 



S9=0 



s\ (20) 



)s 2 



s 2 (21) 

s^=0 



Even though the total energy of the vesicle formally depends on two parameters, 
tq and s 2 , the constraint requiring a constant vesicle volume eliminates one degree 
of freedom, thus yielding ( 1201421 1) . It is also worth noting that the fact that neither 
prolate (s 2 > 0) nor oblate shapes (s 2 < 0) have a bending energy lower than those 
of a sphere means that the expansion for the bending energy d20l contains no linear 
term. 

Equilibrium vesicle deformation, expressed in terms of s 2 , can then be calcu- 
lated by minimizing the total free energy over s 2 : 

j — (Gbend + Gfi e id) = . (22) 

ds 2 

Expressing the equilibrium vesicle deformation s 2 from ( 122b at given condi- 
tions yields: 

' 2 -i^^^^ + TT^ + TT^J • (23) 

As one can see, the dependence of vesicle deformation s 2 on the angular fre- 
quency CO = 2nv contains two dispersion terms. It is worth emphasizing that the 
only approximation used in deriving the expression d23l) is that of a small defor- 
mation (s 2 <C tq). Also worth noting is the fact that the only parameter related to 
membrane elasticity, k c , only scales s%', frequency-dependent electrodeformation 
of vesicles is entirely governed by the electrical parameters of the membrane and 
the aqueous medium and by the vesicle size. 

Fig. [3] shows vesicle deformation in dependence of the angular frequency CO 
of the applied AC electric field. Its most prominent feature is that the deformation 
exhibits a prolate-to-oblate transition at CO ~ 10 4 Hz (i.e., V ~ 10 3 Hz) in the case 
when Oi n < <7 ex (e.g., Oi n = 0.8 a ex , solid curve). On the other hand, vesicle shape 
remains prolate up to CO ~ 10 7 Hz (v ~ 10 6 Hz) when > a ex (e.g., a; n = 1.3 a ex , 
coarsely dashed curve). In both cases, a transition to spherical shape (s 2 ~ 0) is 
observed at CO ~ 10 7 Hz (v ~ 10 6 Hz). All these features are in agreement with the 
experimental findings [181. The transition to spherical shape, however, occurs at as 
low as CO ~ 10 4 Hz when the conductivities inside and outside are equal, — <7 ex 
(Fig-El finely dashed curve). Some of the points in Fig.[TJaround a m / a ex = 1 may 
indicate that such transitions have also been observed. It needs to be pointed out, 
however, that this is not a distinguished property of the a; n = <7 ex case, but merely 
a coincidence at given values of material properties and at given vesicle size. 

The lengthy expressions for the coefficients ^»o, %%, %\ and x 2 figuring in 
( 123b can be somewhat simplified. First of all, it is reasonable to assume that the 
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Fig. 3 Vesicle deformation as a function of the frequency co = 2nv of the applied AC electric 
field. Vesicle deformation is shown both in terms of quadrupolar deformation S2 normalized 
to rg (primary y-axis, left) and vesicle semi-axes ratio c/a ~ (10+^2)/ ( r o — - s 2/2) (secondary 
y-axis, right). The three curves were computed for three different ratios between the electrical 
conductivity of the external aqueous medium and the electrical conductivity of the internal aque- 
ous medium. The dotted line separates prolate shapes from oblate ones. Other parameters used 
for computation were ro = 20 (im, <r ex = 50 (iS/cm, <r m = 10~ 14 S/m, e w = 80£o, e m = 2.5 £0, 
h = 4 nm, k c = 1.2 x 10~ 19 J, E = 500 V/m. 



dielectric permittivity of the aqueous medium inside the vesicle does not differ sig- 
nificantly from the dielectric permittivity of the external medium, £ ex = £ m = £ w 
(this is an approximation; it is known that the dielectric permittivity of water de- 
creases with the increasing concentration of salt in it; e.g., the dielectric permit- 
tivity of physiological saline is only w 72 £0, [20|). A further simplification is 
possible in the case when realistic values for material parameters are taken into 
account: (a) the conductivity of the aqueous solution greatly exceeds that of the 
membrane, ct w S> cr m , and (b) with £ w w 80£o, £ m 2.5£o, and ro/h ~ 1000, then 
£ w *C (ro/h)e m . Retaining only the largest terms, one obtains: 

_ £w(4r c m + £ w ) 
(3r c m + 2£ w ) 2 ' 

(25) 

CT ex 2+X 

r Q c m 2 + x 

T 2 « — — . (26) 

CT ex 2x 

Here, x = Cm/Cex nas been introduced for the ratio of the conductivities of the 
aqueous medium inside and outside the vesicle. 



Frequency-dependent electrodeformation of GUVs 



9 



1.0 



0.0 



-0.5 



0.01 0.1 1 10 100 



Fig. 4 The coefficients £i and £2 figuring in J23b plotted for different values of the ratio of the 
internal and the external electrical conductivity, x = <3i„/o ex . Other parameter values used for 
computation were the same as in Fig. [3] 



The two characteristic times are related to two distinct physical processes: Ti is 
the relaxation time for the Maxwell-Wagner interfacial polarization, and T2 is the 
charging time of a spherical capacitor. In the literature (see, e.g., [21 16,22, 3 1) the 
expressions for them are often encountered in a slightly more general form, which 
distinguishes between the dielectric permittivity of the internal and the external 
aqueous medium. 

As expected, does not depend on conductivities, and the curves for dif- 
ferent Oj n /aex plotted in Fig. [3] merge into a single one at CO — > °°. One can also 
see that in the high-frequency limit, the deformation is always slightly prolate, 
^00 w (h / ro) (e w / £ m ) < 0.01. Such subtle deformation is however in reality likely 
to be subdued by the thermal fluctuations, and is therefore difficult to detect ex- 
perimentally. 

Fig. |4] shows how the coefficients £,\ and ^2 change when the ratio of the in- 
ternal and the external electrical conductivity x = a- m j a ex is varied within two 
decades at given values of material parameters. On the interval x € [i/ioo, 100], £1 
and ^2 can be approximated by rational functions of x: 



(r c m + £ w )(3r c m + 4e w ) + (r c m + £ w )(3r c m + 4e w )x 



L m 



L m 



16(r 2 c2 



L m 



-5r c m e w + 4e2) + (3r ( 
-3r„c m £ w + 24)-(9^c 



,2 2 
C m 



Yhv () c m £ w 
■f- 15r^c m £ w 



I2el)x 
f44)x 



4(3rgc2 + 13r c m e w +44) +4(4^ + 20r o c m e w + 24e2)jc 



(27) 
(28) 



It is possible to derive an approximate expression for the critical value of x, 
below which the prolate-to-oblate transition occurs. In the intermediate region, 
1/t 2 < CO < 1/ti, the function %{oo;x) = + £1/(1 + co 2 z 2 ) + 1 2 /(1 + co 2 x\) 
figuring in ( 123b approximately evaluates to % {C0\x) m ^ M + £1. This needs to be 
negative in order for an oblate deformation. From this requirement, an approxi- 
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Fig. 5 Phase-contrast micrograph of a vesicle in an AC electric field in response to the frequency 
of the applied field: (a) 1 kHz, (b) 1 1 kHz, (c) 29 kHz. The direction of the electric field corre- 
sponds to the horizontal direction. The bar represents 10 urn. The conductivity of the aqueous 
medium is Ojn =17 (iS/cm. 



mate expression for jc cr j t can be obtained: 

27r c m + 93e w 



27r c m + llle„ 



(29) 



Substituting realistic values for material parameters and vesicle size into ( I29K one 
finds that x cr ; t sa 0.97 for a vesicle with r$ = 5 (im, x cr j t « 0.986 for a vesicle with 
ro = 10 um, and x cr i t « 0.997 for a vesicle with ro = 50 urn. 



3 Experiment 

3.1 Vesicle preparation 

Giant unilamellar vesicles (GUVs) were prepared from commercially available 
(Avanti Polar Lipids, Alabaster, AL, USA) synthetic 1 -palmitoyl-2-oleoyl-.m-gly- 
cero-3-phosphocholine (POPC) in double-distilled water using the electroforma- 
tion method [23 24 1 in a chamber which allows for an easy access of the resulting 
vesicle suspension [25 1. In order to alter the conductivity of the solution, sodium 
chloride was added to the solution in concentration of up to 0.2 x 10~ 3 mol/L. 



3.2 Experimental procedure 

The experimental chamber consisted of two 0. 1 mm thick stainless steel electrodes 
mounted on an object glass, leaving a 0.6 mm wide gap between the electrodes. 
The cover slip was mounted with vacuum grease (Baysilone; Bayer, Leverkusen, 
Germany), enclosing a few drops of vesicle suspension. The chamber was placed 
onto an inverted light microscope (Zeiss/Opton IM 35, objective Zeiss Ph2 Plan 
40/0.60; Zeiss, Oberkochen, Germany), and phase contrast technique was used. 
The micrograph was recorded using a CCD camera (Cohu 6700; Cohu, San Diego, 
USA), taped using a U-Matic VCR (Sony VO-9800P) and later digitized on a PC 
with a frame grabber (Matrox Meteor II; Matrox, Dorval, Canada). The same com- 
puter was used to set the voltage and frequency of the applied electric field, driving 
a GPIB -controlled function generator (Iskra MA-3735; Iskra, Horjul, Slovenia). 
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Fig. 6 Prolate-to-oblate transition of a phospholipid vesicle in the kHz range. Vesicle deforma- 
tion is expressed as the ratio of its semi-axis in the direction along the field (c) and its semi-axis 
perpendicular to the field (a). Vesicle radius rg = 19.8 ±0.1 urn, conductivity of the aqueous 
medium <7j n = 4.6 (iS/cm. 



Using an acoustic coupler connected to RS-233 port and a telephone handset with 
an audio jack, the data about the selected voltage and frequency were simultane- 
ously recorded onto one audio channel, while the other audio channel was left for 
the experimentalist's comments during the course of experiment [26|. 

In the experiment, the prolate-to-oblate transition frequency for a given vesicle 
was first coarsely estimated by testing the effect of the electric field at 2-3 differ- 
ent frequencies. After this coarse estimate, a programmed sequence of step-wise 
frequency change was applied. The frequency step ranged from 100 Hz to 1 kHz, 
depending on the vesicle, and the duration of the step was approximately 3 s. This 
duration was experimentally selected as being long enough for the vesicle to adapt 
to a small frequency change, yet short enough to minimize dielectrophoretic drift 
and allow repetitive measurements. The programmed sequence was repeated sev- 
eral times with increasing and decreasing frequencies in order to observe possible 
hysteresis (Fig. [6]). To allow a more precise determination of transition frequency, 
the range was narrowed and the frequency steps were decreased in the later runs 
(symbols □, A, and y in Fig.[6]l. 



4 Results 

4. 1 Analysis of the experimental data 

A total of 46 vesicle recordings were examined. After discarding the recordings 
with a vesicle relative volume too close to 1, where no noticeable shape change 
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Fig. 7 The dependence of the prolate-to-oblate transition frequency on the vesicle size and the 
conductivity of the aqueous medium. Three runs of experiments with three different conduc- 
tivities of the aqueous medium were performed: 17, 4.6, and 1.3 uS/cm. The solid line corre- 
sponds to numerically solving si = (Eq,l23t for to, with a m = 17 um, and the dashed line with 
(Ti n = 4.6 (im. In both cases, t7i n /(T ex = 0.9. Other parameter values used for computation were 
the same as in Fig. [3] 



was observed, as well as the recordings in which the experiment was interrupted 
by the dielectrophoretic drift of a vesicle to a region where further observations 
were not possible, 21 measurements on different vesicles were taken into con- 
sideration for analysis. Of these, 8 were prepared in a medium with conductivity 
17 uS/cm, 3 in a medium with conductivity 4.6 uS/cm and 10 in a medium with 
conductivity 1.3 uS/cm. The transition frequency for the prolate-to-oblate transi- 
tion vs. vesicle size is plotted in Fig. [7] Transition frequencies are reciprocally 
related to the vesicle radius, which is consistent with the recently published find- 
ings of another group (K. Antonova et ah, in press). 

The existence of prolate-to-oblate transitions indicates that the ratio Oi n /<7 ex 
in the experiments was below x cr it, even though the compositions of the aqueous 
solution in the vesicle interior and the vesicle exterior were initially identical. We 
would like to propose an explanation for this phenomenon. First, we need to em- 
phasize that the conductometer electrode we used requires at least 6 ml of sample, 
therefore it was impossible to measure directly the conductivity of the vesicle so- 
lution. Instead, the conductivity of the aqueous medium was measured before it 
was used for hydrating the phospholipid film. The conductivity of the aqueous 
solution may increase afterwards by the impurities introduced either at the elec- 
troformation process, during storage or during sample preparation. In order to test 
this hypothesis, we performed a separate experiment, where we simulated the ma- 
nipulation of a vesicle sample. Two flasks were filled with 0.03 mmol/L solution 
of NaCl in 0.1 mol/L sucrose solution (a = 14.5 uS/cm). The first flask was kept 
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sealed in cold storage during the course of the experiment, while the second flask 
was brought to room temperature every day, and a small amount of the solution 
was pipetted out for conductivity measurement before returning the flask back to 
cold storage. While the conductivity of the solution in the first flask remained un- 
changed, the conductivity of the solution in the second flask increased by w 8% 
with each consecutive pipetting, reaching 31.5 (iS/cm after 11 days. We can pre- 
sume that similar processes occur in vesicle suspension, where the conductivity of 
the solution in the vesicle interior remains unchanged as the phospholipid mem- 
brane is virtually impermeable for ions, while the conductivity of the solution in 
the vesicle exterior increases with each successive pipetting. Since even for the 
smallest of the vesicles shown in Fig.[7](ro sa 5 um), Xcdt w 0-97 (cf. Eq.[29]l, it is 
clear that at such low conductivities of the solution, an 8% increase of the conduc- 
tivity of the external medium caused by a single pipetting is sufficient to bring the 
vesicle suspension into the regime where oblate shapes exist in the intermediate 
frequency range. 

The two computed lines in Fig. [7] were obtained by numerical root-finding of 
S2 = 0, S2 being defined by ( 123b . The initial conductivity of the solution (17 uS/cm 
and 4.6 uS/cm, respectively) was taken as the conductivity of the internal aqueous 
medium, while the conductivity of the external aqueous medium is 1/0.9 a; n , or 
approximately 10% higher. The quantitative agreement with the measured data 
at 17 uS/cm and 4.6 (iS/cm is within the experimental error, while the data for 
1.3 (iS/cm coincide with those for 4.6 (iS/cm. This may indicate that the actual 
electrical conductivity Oj n in this latter case was not 1.3 (iS/cm, but somewhat 
higher (estimated around 5 uS/cm). 

4.2 Morphological diagram 

In order to plot a morphological phase diagram akin to the one shown in Fig. Q] 
two distinct regimes have to be considered. For <7i n /a ex < 1, vesicle deformation 
$2 ( 1231) changes sign, and the prolate-to-oblate and the oblate-to-spherical transi- 
tions can be computed by solving s% = numerically. For (7; n / a ex > 1, S2 remains 
positive, and instead one has to rely on the characteristic times ( 125 1 126| |. Here, 
the relaxation term containing Ti corresponds to prolate-to-spherical transition in 
Fig.[T] The relaxation term corresponding to T2 corresponds to an unobserved tran- 
sition, i.e., a prolate vesicle is expected to become slightly less prolate (see also 
the dashed curve in Fig. O at frequencies CO > I/T2. Such subtle changes, how- 
ever, are probably difficult to quantify experimentally, and have therefore not been 
reported in [18 |. 

The results, shown in Fig.[8l show a reasonable agreement with the experimen- 
tal diagram (Fig. [TJ. In particular, for a; n < a ex , the prolate-to-oblate transition 
frequency increases with increasing <7i n /a ex . The prolate-to-spherical transition 
for Oi n > <7 e x a l so increases with increasing Oi n /(7 ex , while the oblate-to-spherical 
transition, observed for (J m < a ex , decreases with increasing a- m / a ex . All these 
observations are in agreement with the experiment. 

We would also like to comment on the reported coexistence of prolate and 
oblate vesicle shapes at 1 w 1 in the intermediate frequency range II 1 81 . As we 
have shown earlier (Fig. 13, the prolate-to-oblate transition frequency exhibits 
a strong dependence on vesicle size: the transition frequency for a fairly large 
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to [Hz] 

Fig. 8 Prolate and oblate vesicle deformation plotted for different values of the ratio of the 
internal and the external electrical conductivity, x = a m /a ex . In the x < 1 region, the sign changes 
in <23> . and either prolate-to-oblate transition or oblate-to-spherical transitions are computed, 
while forx > 1, reciprocal characteristic times 1/Ti and I/T2 are plotted. Other parameter values 
used in evaluation were the same as in Fig. [3] The dashed curve corresponds to the relaxation 
which does not manifest itself as a morphological shape transition. 



vesicle (ro = 50 urn) is V2 ~ 3 kHz, while for a small vesicle (ro = 5 urn) it is 
V2 ~ 30 kHz. Within this frequency interval, it is thus quite feasible to observe a 
small prolate vesicle simultaneously with a large oblate vesicle. While the authors 
do not comment on the relative sizes of the vesicles in question, we believe this 
phenomenon can be explained if: (a) the a; n / a ex ratio was below x cr i t (i.e., slightly 
below 1), thus making oblate shapes possible in the intermediate frequency range, 
and (b) the frequency was above the prolate-to-oblate transition frequency for the 
larger vesicle and below the prolate-to-oblate transition frequency for the smaller 
vesicle. 



5 Discussion 

5.1 Comparison with earlier results 

While the model presented above builds on the legacy of the model developed 
over 20 years ago by Winterhalter and Helfrich 1121 . it departs from it in two 
aspects. Firstly, it allows for a different electric conductivity in the vesicle inte- 
rior and the vesicle surroundings, and secondly, it is simplified by treating the 
vesicle membrane as a thin shell right from the start. The first modification was 
a necessity imposed by the experiments presented in [18|. The limitation to thin 
shells is appropriate since any results applicable to giant vesicles fall into this 
limit (7 sa 1 in 1121 "). The results presented here thus do not differ significantly 
from those obtained by the authors who started from a general finite-thickness 
model (T. Yamamoto et ah, to appear). Treating membrane thickness as finite is, 
however, always required when when one attempts to tackle membrane dielectric 
anisotropy EEH^EDED- 
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Hyuga et al. ifTfjl provide treatment both for the vesicle deformation in an 
AC electric field and the response of the vesicle to a step-wise change in E, and 
arrive at the expressions (1251126b . However, they later conclude that the fast mode 
dampens very quickly and that only the effects of the slow mode are expected to 
show a visible influence on deformation dynamics, and consequently focus on the 
lower (i.e., prolate-to-oblate) transition. 

More recently, a theoretical analysis of the experimental results presented in 
[ 18 1 has been published by the same group 1321 . where a force-balance approach 
is used instead of energy minimization. While this does not introduce any major 
advantage in the problems treated here, i.e., the phase diagram of vesicle shapes, 
it allowed a more correct prediction of not only the nature of deformation, but also 
of its extent (see below). Even more important, by extending the treatment beyond 
the equilibrium shape it opens a possibility of treating vesicle dynamics. 

Finally, we find it necessary to relate the results in this paper with our own 
earlier results [31], where we argued that the prolate-to-oblate transition can be 
explained by the dielectric anisotropy of the membrane. In view of later exper- 
imental findings III 81 . which put the prolate-to-oblate transition into a broader 
context, as well as of the fact that even a very minute difference in the electrical 
conductivities of the aqueous medium inside and outside the vesicle (cf. eq. |29| l 
makes prolate-to-oblate transition possible, we believe it is quite possible that the 
observed prolate-to-oblate transition was caused primarily by a minute unintended 
increase of the conductivity of the external aqueous medium, which may have 
overshadowed other possible effects. Even though the frequency-dependent elec- 
trodeformation of giant vesicles is likely to be primarily governed by other effects, 
we however continue to believe that the dielectric anisotropy of the phospholipid 
membrane is a phenomenon worth investigating. 

5.2 Quantitative estimate of vesicle deformation 

Though simple, the model presented in this paper provides a fairly accurate de- 
scription of the experimentally obtained morphological phase diagram (Fig. [TJ, 
i.e., it predicts the nature of deformation (prolate, oblate, or spherical), and the 
frequencies at which the transitions between these shapes occur. However, being 
limited by design to deformations close to spherical, this model, like the original 
Winterhalter-Helfrich model, does not predict the extent of deformation correctly. 
In the case of the Winterhalter-Helfrich model, this feature has already been noted 
01811321 . For the most part, the extent of deformation is determined by the relative 
volume of the vesicle, i.e., a flaccid vesicle deforms more than a nearly spheri- 
cal one. For simplicity, the treatment presented in this paper omits the effects of 
thermal fluctuations, which have already been treated elsewhere 1133 11341 . 

5.3 Application in microliter conductometry 

A quantitative model of vesicle shape with respect to the electrical conductivity of 
the medium inside and outside the vesicle allows us to exploit this phenomenon 
from the other end: a vesicle with a known internal conductivity and a measurable 
size can serve as a probe for sensing the electrical conductivity of the medium 
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in its surroundings. For this potential application, the kHz-range prolate-to-oblate 
transition seems particularly well-suited, because (a) it is an easily observable 
qualitative change in the vesicle shape which occurs in a narrow frequency in- 
terval, and (b) the transition frequency exhibits a strong dependence on vesicle 
size, which allows for some tuning, i.e., in a heterogeneous vesicle population, 
one can select the vesicle which manifests the effect most prominently. Recent 
developments in vesicle preparation methods [35 36] also bring the physiological 
conductivities into reach. 



6 Conclusions 

It has been demonstrated that a modified Winterhalter-Helfrich model lfl2l which 
obtains a parametrized vesicle shape by minimizing the sum of the membrane 
bending energy and the energy due to the electric field, can be successfully ap- 
plied to (a) quantitatively explaining the experimentally observed prolate-to-oblate 
shape transition in the kHz region, and (b) more generally, to the experimentally 
obtained morphological phase diagram of GUVs with respect to two parameters, 
the frequency of the applied AC electric field, and the ratio of the electrical con- 
ductivities of the aqueous media inside and outside the vesicle [ 18 1. The obtained 
results are compared with the findings in the literature, and the limitations of the 
model concerning the quantitative prediction of membrane deformation are dis- 
cussed. Finally, an application of the observed effects in the conductometry of 
microliter samples is proposed. 
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